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Abstract 


This  paper  describes  an  analysis  of  6.4  to  25.4  mm 
(0.25  to  1.0  in.)  thick  composite  laminates  subjected  to  uniaxial 
compressive  loading  that  experimentally  showed  a  decrease  in 
strength  with  increasing  thickness.  The  analysis  was 
performed  to  determine  if  the  reduction  in  stretch  was  an 
intrinsic  thickness  material  effect  or  if  it  could  be  attributed  to 
through-thickness  restraint  on  the  specimen  caused  by  the  test 
fixture.  The  analysis  was  based  on  closed  form  solutions  for  the 
formation  of  kink-band  failures  in  the  presence  of  fibers 
misaligned  with  the  principal  axis  of  compression  loading. 

The  fiber  misalignment  was  determined  by  a  finite  element 
analysis  that  accounted  for  the  displacement  of  the  laminate 
outer  plies  where  the  laminate  exited  the  compression  test 
fixture.  The  correlation  between  the  experimental  results  and 
the  theoretical  analysis  showed  the  compression  strength  of  the 
AS4 1 3501-6  and  S2 13501-6  laminates  to  he  independent  of 
thickness  and  directly  proportional  to  the  through-thickness 
fixture  restraint  on  through-thickness  Poisson  expansion.  This 
theoretical  and  experimental  comparison  also  demonstrated  a 
strong  correlation  between  fiber  misalingment,  its  resulting 
shear  stress  state  and  the  kink-band  compression  failure 
mechanism. 
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Introduction 

Many  factors  influence  the  compressive  response  of  composite 
materials  and  considered  together  or  separately  they  can  trigger  a  number 
of  failure  modes  (1).  One  factor  that  has  been  shown  theoretically  to  have  a 
significant  effect  on  the  compression  strength  of  composite  materials  is 
fiber  misalignment  or  fiber  curvatizre  (2-6).  Either  of  these  conditions 
result  in  reinforcing  fibers  that  are  at  an  angle  to  the  axis  of  primary 
compressive  load.  Shear  stresses  can  become  very  high  even  at  small 
misalignment  angles.  In  addition,  carbon  and  fiberglass  reinforced 
composites  exhibit  nonlinear  shear  stress-strain  response,  resulting  in  a 
diminishing  shear  moduli  with  increasing  shear  deformation.  Coupled 
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with  the  observation  that  microbuckling  based  compression  failure  theories 
show  a  strong  relationship  between  shear  modulus  and  compression 
strength  (6,  7),  it  is  imderstandable  that  small  fiber  mislignments  could 
significantly  affect  compression  response  by  introducing  moderate  to  high 
shear  stresses  in  the  material.  Although  microbuckling  failure  is  not  the 
only  compression  failure  mechanism  possible,  and  has  not  been 
conclusively  shown  to  be  the  most  predominant  failure  compression 
mechanism  in  composites,  many  investigators  have  experimentally  shown 
a  direct  relationship  between  matrix  composite  shear  modulus  and 
iiltimate  compression  strength  (8,  9). 

One  major  weakness  of  the  hypothesis  that  compression  response  is 
strongly  influenced  by  the  presence  of  fiber  misalignment  and  curvature  is 
the  lack  of  experimental  data  showing  direct  correlation  between  actual 
fiber  misalii.’nment  or  curvature  and  compression  failure  initiation.  The 
observation  that  a  fiber  misalignment  of  1®  to  3®  from  principal  load 
direction  significantly  reduces  compression  theroretically,  and  the 
assertion  that  there  are  fibers  misalingned  to  this  degree  in  every  composite 
compression  specimen  does  not  provide  substantiation  to  the  subject 
hypothesis. 

In  this  paper  a  correlation  between  experimentally  determined 
compression  strength,  fiber  misalignment  angle,  and  theoretical 
compression  strength  is  presented.  This  paper  describes  the  details  of  the 
theoretical  portion  of  this  correlation  with  reference  to  the  experimental 
portion  that  deals  with  the  compression  response  of  thick  (6.4  -  25.4  mm 
[0.25  -  1.0  in.])  laminates  (10).  A  schematic  diagram  of  the  fixture  used  for 
the  compression  testing  is  shown  in  Figure  1.  A  principal  objective  of  the 
experimental  program  was  to  determine  the  mechanical  characteristics  of 
thick  composites  in  compression  and  observe  if  they  differ  from  those  seen 
in  thinner  composite  materials.  The  conclusions  drawn  from  this 
experimental  program  are: 

-  the  compressive  elastic  constants  for  AS4/3501-6  and  S2  glass/350 1- 
6  do  not  change  as  specimen  thickness  increases, 

-  however,  significant  through-thickness  strains  due  to  nonlinearities 
in  Poisson's  ratio  developed. 
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-  there  was  a  20%  drop  in  [02/90]ns  compression  strength  when 
specimen  thickness  increased  from  6.4  to  25.4  mm  (0.25  to  1.0  in.), 

-  the  failure  mechanisms  did  not  change  with  increasing  thickness 
and  were  predominantly  kink  band  formation  and  delamination,  as 
seen  in  thin  composite  materials, 

-  compression  failures  occurred  predominantly  at  the  gage-section/ 
tab  termination  interface. 

-  in  autoclave  cxared,  high  quality  AS4/3501-6  and  S2  glass/3501-6 
[02/90]ns  laminates,  through-thickness  lamina  level  waviness  along 
the  0°  axis  was  0.660®  and  0.820°  respectively,  ref.  (11). 

Figure  2  shows  kink  bands  at  the  tab  termination  region  of  a  48  ply,  6.4  mm 
(0.25  in.)  [02/90]8s  AS4/3501-6  specimen  after  compression  failure. 

The  fact  that  the  strength  of  the  AS4  and  S2  glass  laminates  decreased 
20%  even  though  all  other  characteristics  remained  unchanged  suggest 
this  phenomenon  should  be  further  investigated  to  determine  if  it  is 
primarily  due  to  material  scaling  effects  or  fixture  imposed  effects.  Since 
most  of  the  laminate  compression  failures  occurred  in  the  region  where  the 
tabs  terminated,  the  effect  of  experimental  technique  on  failure  was 
theoretically  investigated  and  is  presented  in  this  paper.  In  this  region  of 
the  specimen,  fiber  curvature  due  to  the  combined  effect  of  gage  section 
Poisson  expansion  and  through-thickness  fixture  restraint  was  the  largest 
(Figure  3). 

In  consideration  of  the  above  observations,  the  objective  of  the  present 
analysis  was  to  determine  if  the  reduction  in  compression  strength 
measured  in  the  series  of  tests  on  thick  coupons  could  be  attributed  to  the 
total  through-thickness  fiber  misalignment.  Initial  fiber  waviness  was 
directly  measured  while  misalignment  due  to  fixture  restraint  was 
predicted  through  finite  element  analysis.  This  information  was  then 
substituted  into  existing  expressions  for  compression  strength  modified  to 
include  both  components  of  fiber  misaligmnent  under  consideration  and  to 
include  stress  concentration  effects  from  the  fixture  induced  through- 
th’ckness  restraint.  While  the  analysis  reflects  specific  conditions  defined 
b  le  experimental  program,  this  study  points  out  the  increasing  influence 
of  nxture  restraint  on  inplane  compression  strength  as  section  thickness 
increases. 
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Theoretical  Procedure 

Since  compression  strength  has  been  shown  theoretically  to  be 
strongly  influenced  by  fiber  curvature  and  misalignment,  and  the  value  of 
Vxz  measured  for  the  thick  laminates  is  large,  the  amount  of  effective 
thickness  expansion  (z-direction)  in  the  gage  section  was  examined.  The 
effective  through-thickness  gage  section  expansion  is  defined  as  the 
difference  between  the  z-direction  expansion  that  occurs  within  the  gage 
section  of  the  specimen  and  that  which  occurs  within  the  clamped  tabbed 
region  of  the  specimen  during  the  application  of  a  compressive  load.  The 
effective  gage  section  expansion  is  shown  schematically  in  Figure  3.  This 
expansion  is  of  concern  since  it  results  in  a  fiber  curvature  where  the 
coupon  exits  the  fixture  clamping  blocks,  and  is  maximiim  on  the 
specimen  outer  plies.  Figure  3  also  shows  the  angle  at  which  the  specimen 
outer  ply  exits  the  clamping  blocks  as  a  result  of  this  expansion.  If  the 
effective  gage  section  expansion  increases  with  specimen  thickness,  then  so 
will  fiber  curvature,  resulting  in  theoretically  lower  compression  strengths 
for  thicker  composites. 

The  three  displacement  components  that  influence  effective  gage 
section  expansion  are  the  free  gage  section  expansion,  the  specimen 
contraction  under  the  clamping  blocks  due  to  bolt  preload,  and  the 
specimen  and  tab  expansion  within  the  fixture  clamping  blocks,  and  are 
shown  in  Figure  3.  In  order  to  quantify  these  effects,  the  following 
theoretical  approach  was  applied: 

1)  determine  effective  gage  section  expansion  tising  the  strength-of- 
materials  bolt  analysis  and  general  Hooke's  Law  for  an  orthotropic 
material, 

2)  use  this  expansion  as  a  uniform  displacement  boundary  condition  in 
a  finite  element  analysis  to  determine  outer  ply  displacement 
geometry, 

3)  determine  the  effect  of  fiber  misalignment  on  compression  strength 
using  a  compression  failure  theory  that  is  based  on  a  kink  band 
failure  mechanism. 
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In  addition  to  this  approach,  a  finite  element  model  of  1/4  of  the 
complete  fixture-specimen  assembly  with  uniform  end  loading  was  also 
performed.  The  gage  section  expansion  determined  using  this  approach 
was  the  same  as  that  determined  in  step  1)  above.  Since  a  large  number  of 
elements  were  needed  for  the  convergence  study  performed,  the  strength-of- 
materials  solution  for  gage  section  expansion  was  utilized  and  only  a 
portion  of  the  specimen/tab  for  each  laminate  thickness  region  was 
analyzed  using  the  finite  element  procedure. 

The  free  gage  section  expansion  can  be  expressed  as: 


CFcUVxz 


where 

Gc  =  longitudinal  compressive  stress  at  specimen  failure 
ts  =  specimen  thickness  in  the  gage  section 
Vx2  =  through-thickness  Poisson's  ratio 
Ex  s  longitudinal  modulus 


Preload  contraction  can  be  expressed  as: 


r^t 


where 
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Oz  =  through-thickness  compressive  stress  due  to  bolt  preload 
t  =  specimen  thickness  within  the  tabbed  section  (including  tabs) 

Ez  =  through  thickness  modulus 

and  the  specimen  and  tab  expansion  within  the  fixture  clamping  blocks  due 
to  the  applied  longitudinal  compressive  stress  as: 
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where 

Ob  =  bolt  stress  at  specimen  failure,  reference  (12) 
Lb=  bolt  length 
Eb  =  bolt  modiolus 


Finally  the  effective  gage  section  expansion,  Agfr  can  be  expressed  as: 


.  actsVx7  ^ 
^etf=  — t -  + 
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The  finite  element  code  used  to  determine  the  outer  ply  displacement 
geometry  resulting  from  the  gage  section  expansion  was  ABAQUS 
developed  by  HKS,  and  was  nm  on  a  Digital  Equipment  Corp.  MicroVax.  A 
2-D  analysis  using  2nd  order  plane-strain  orthotropic  elements  for  the 
composite  specimen  and  tabs  was  performed. 

From  the  finite  element  modeling  the  outer  ply  displacement  and  exit 
angle  (shown  in  Figure  3)  were  determined.  This  information  was  then 
used  as  input  into  failure  theories  that  are  based  on  kink  band  formation  as 
the  principal  mechanism  of  failure  for  fiber-reinforced  composite 
materials. 

Argon  (2)  8ind  Budiansky  (3)  have  both  proposed  compression  failure 
theories  that  were  developed  based  on  the  observation  of  kink  band 
formation  in  compressively  loaded  orthotropic  materials.  They  also  include 
a  fiber  misalignment  angle  term  in  their  analysis  as  the  geometric 
perturbation  that  creates  a  region  of  increased  shear  stress  and  triggers 
failure  due  to  instability. 

In  the  work  by  Argon,  he  considers  a  region  of  initial  misalignment  ( 
^o)  with  respect  to  the  loading  axis  in  a  fiber  reinforced  composite.  This 
region  causes  an  interlaminar  shear  stress  and  a  corresponding  element 
rotation,  that  further  increases  shear  stress  until  shear  collapse  instability 
occurs.  This  model  parallels  an  analysis  of  the  formation  of  kink  bands  in 
metal  crystals.  His  derivation  leads  to  an  expression  for  the  compression 
strength  of 
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(X^)a  =  ^ 
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where  is  the  compression  strength,  S13  is  the  ultimate  transverse  shear 
strength,  and  ((>0  is  the  initial  misalignment  angle. 

Budiansky  extended  Rosen's  (7)  elastic  shear  microbuckling  solution 
for  the  ultimate  compression  strength 

(Xi)b=Gi3  6 

to  include  plasticity,  assiiming  perfect  plasticity  in  shear  beyond  y=  E13  = 
S13/G13  and  fiber  misalignment.  The  resulting  expression  is 


The  analj  ses  by  both  Argon  and  Budiansky  assume  the  presence  of 
an  initial  fiber  misalignment.  In  applying  these  analyses  to  the  thick 
compression  specimens  tested  in  this  program,  fixture  induced  fiber 
misalignment  that  occurs  with  the  application  of  a  compressive  load  is 
considered  equally  as  important  as  the  initial  misalignment.  Therefore  the 
term  (|)o  in  v  qs.  5  and  7  is  generalized  to  include  both  initial  misalignment 
(0w  since  in  this  study  lamina  waviness  creates  the  misalignment)  and 
Poisson  induced  misalignment  By.  In  the  present  analysis  the  magnitude 
of  oi  in  the  0°  plies  of  the  [02/90]n8  laminate  is  determined  from  classical 
lamination  theory. 

These  analyses  eJso  do  not  include  the  effects  of  transverse  inplane 
and  through-thickness  normal  stresses.  For  the  [02/90]ns  laminates  in  this 
study  inplane  transverse  normal  stresses  are  present  from  transverse  ply 
restraint  and  through-thickness  normal  stress  from  fixture  restraint. 
Inplane  normal  stresses  are  neglected,  and  through-thickness  normal 
stress  contributions  are  accoimted  for  with  the  addition  of  a  stress 
concentration  term  on  the  longitudinal  lamina  stress  that  is  determined 
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from  laminate  plate  theory  and  far  field  laminate  stress.  The  resiilting 
expressions  for  compression  strength  are: 


k  -  ^13 

(««+  0v) 


(8) 


k,(x')b  = 


'13 


.(e.+  ej 


+  e. 


'13 


(9) 


This  stress  concentration  factor  was  defined  through  a  correlation  between 
the  experimental  results  at  a  single  specimen  thickness  (6.4  mm)  and  the 
finite  element  stress  analysis  for  this  thickness.  Consequently  eqs.  8  and  9 
are  not  newly  derived  expressions  for  compression  strength,  but  rather  are 
variations  of  the  original  solutions  by  Argon  and  Budiansky  to  accotmt  for 
observations  from  the  experimental  thick-section  compression  test  results. 


Results  and  Discussion 

The  effective  gage  section  expansion  determined  by  equation  4  is  a 
linear  function  of  specimen  thickness.  The  contributions  from  each  term 
in  equation  4  along  with  the  effective  gage  section  expansion  are  listed  in 
Table  1  for  each  material  and  thickness. 

The  laminate  material  properties  used  in  the  finite  element  analyses 
and  the  effective  gage  section  expansion  calculations  are  the  values  from 
Table  2.  The  theoretical  values  were  determined  using  experimentally 
determined  lamina  properties  (Table  3)  and  Pagano's  exact  theory  for 
determining  3D  orthotropic  laminate  properties  (13,  14).  The  details  of  the 
theoretical  procedure  are  included  in  ref.  (12)  The  values  for  laminate 
shear  moduli  are  initial  tangent  moduli  values,  however  preliminary  finite 
element  analyses  with  shear  moduli  decreased  40%  showed  a  20%  increase 
in  outer  ply  exit  angle.  Therefore  the  shear  moduli  used  in  the  ABAQUS 
finite  element  runs  were  60%  of  those  in  Table  2.  A  value  of  60%  was 
chosen  since  this  moduli  is  equal  to  the  secant  modulus  defined  from  0 
strain  up  to  the  point  where  the  curve  flattens  out  for  the  two  materials 
tested  by  the  ASTM  D3518  method  .  Data  documenting  the  specific  effect  of 
reduced  shear  moduli  on  outer  ply  exit  angle  will  follow. 


8 


In  order  to  provide  the  outer  ply  exit  angle  needed  for  the  models  by 
Argon  and  Budiansky,  the  finite  element  modeling  scheme  shown  in 
Figures  4  and  5  was  used.  These  figures  show  the  finite  element  mesh 
discretization  used  for  the  coarsest  mesh  studied  and  the  application  of  a 
constant  displacement  boundary  condition,  respectively.  The  magnitude  of 
the  through-thickness  displacement  to  be  applied  was  determined  by  taking 
one  half  of  the  effective  expansion  in  Table  2  and  multiplying  it  by  the  ratio 
of  the  free  gage  section  thickness  plus  tab  thickness  to  the  free  gage  section 
thickness.  2nd  order,  plane  strain  elements  with  an  aspect  ratio  of  2:1  were 
used  for  the  analysis.  Although  the  figures  in  this  paper  show  a  line  of 
distinction  between  the  tabs  and  the  specimen,  the  tabs  used  in  the 
experimental  study  were  made  of  the  same  material  and  stacking  sequence 
as  the  specimens  so  all  analyses  assumed  continuous  material  properties 
across  the  tab-specimen  botmdary. 

Figure  6  shows  preliminary  results  for  the  outer  ply  displacement 
geometry  for  the  all  three  thicknesses  of  the  AS4/3501-6  specimens,  and 
shows  an  increase  in  the  outer  ply  exit  angle  with  increasing  specimen 
thickness.  After  these  preliminary  results,  a  mesh  refinement  study 
beginning  with  the  mesh  discretization  shown  in  Figure  5  was  conducted. 
Analyses  with  mesh  refinement  to  1/16  the  size  shown  in  Figure  6  were 
conducted  to  determine  solution  convergence.  The  number  of  elements  and 
corresponding  number  of  nodes  used  in  the  convergence  study  is 
summarized  in  Table  4. 

Table  5  lists  the  values  of  6v  determined  for  each  material  and 
specimen  thickness  and  for  each  mesh  size  analyzed.  A  mathematical 
analysis  of  this  data  indicates  9v  does  not  converge  with  decreasing  mesh 
size  (12).  Consequently,  and  average  exit  angle  is  defined  consistent  with 
the  mesh  discretization  and  the  heterogeneity  of  the  laminated  system  (ie.  a 
repeating  sublaminate).  A  convergent  exponential  form  for  By  as  a 
function  of  z  was  formulated.  This  form  is 


By  =  kz-m  (10) 

The  coefficients  of  this  equation  were  determined  firom  a  plot  of  ln(0n) 
versus  ln(z).  The  values  for  k  and  m  are  listed  in  Table  6.  Integrating 
equation  10  will  yield 
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(11) 


which  is  an  average  value  of  0v  over  a  specified  distance  1. 

The  6v  data  to  be  used  in  the  averaging  procedure  was  taken  from  the 
ABAQUS  runs  with  an  x-direction  element  size  of  0.2  mm  (0.0078  in.)  The 
values  for  the  outer  ply  exit  angles  calculated  using  k  and  m  from  Table  6 
and  equation  11  are  hsted  in  Table  7.  The  dimension  G)  over  which  the 
values  of  6v  were  determined  was  0.38  mm  (0.015  in.)  This  dimension  is 
eqmvalent  to  three  ply  thicknesses,  and  equals  the  dimension  of  the 
repeating  [02/90]  sublaminate  used  in  the  thick  laminates  tested.  This  table 
also  lists  the  average  lamina  level  waviness  measured  for  the  [02/90]ns 
laminates,  and  the  total  fiber  misalignment  angle  resulting  from  a 
superposition  of  the  initial  misalignment  and  the  misalignment  at  failure 
due  to  fixturirg  effects. 

As  previously  mentioned,  the  shear  modtili  used  in  the  calculation  of 
the  outer  ply  exit  angles  in  Table  7  were  60%  of  those  listed  in  Table  2. 

These  values  were  used  since  0v  was  foimd  to  be  sensitive  to  the  laminate 
shear  modudi.  This  sensitivity  was  quantified  using  the  averaging 
technique  for  0.38  mm  (0.015  in.)  of  the  laminate  thickness  as  described 
above  for  the  calctdation  of  the  outer  ply  exit  angle  with  laminate  thickness. 
The  effect  of  shear  moduli  nonlinearity  are  shown  in  Table  8,  and  the 
results  show  a  25%  increase  in  outer  ply  exit  angle  with  a  40%  drop  in 
shear  moduli.  The  actued  effect  of  shear  nonlinearity  may  be  greater  than 
shown  here,  since  only  a  linear  analysis  was  performed,  and  the 
instantaneous  shear  moduli  that  could  be  used  in  a  nonlinear  analysis 
would  be  much  less  them  the  secant  moduli  used  here.  A  nonlinear 
analysis  could  be  performed  by  exercising  the  UMAT  option  in  ABAQUS 
that  allows  for  nonlinear  material  properties  through  user  written 
material  property  subroutines.  Although  not  used  for  this  program,  this 
option  is  being  investigated  in  related  Navy  programs. 

Considering  equations  8  and  9,  (6w  +  By)  is  now  known  and  has 

been  determined  experimentally.  The  laminate  strength  for  each  specimen 
thickness  are  listed  in  Table  9.  Since  the  stress  at  failure  in  the  outer  0° 
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plies  is  necessaiy  for  these  equations,  the  0®  ply  Xj  (Table  9)  was 

determined  using  classical  lamination  theory  and  the  lamina  properties 
from  Table  3,  including  residual  thermal  stresses  for  a  AT  of  138®C. 

Si3  and  ei3  are  also  required  to  evaluate  kt  from  equations  8  and  9. 
From  the  [±45]  tension  tests  run  to  determine  G13  in  Table  3,  S13  was  fovmd 
to  be  75.8  MPa  (11,000  psi)  for  both  the  AS4/3501-6  and  the  S2  glass/3501-6 
materials.  Coupon  failure  did  not  occur  at  75.8  MPa  (11,000  psi)  for  these 
materials,  but  the  shear  stress-strain  curve  flattened  out  at  this  stress.  To 
be  consistent  with  Budiansky's  (7)  assximption  of  elastic  perfectly  plastic 
shear  response  beyond  a  given  shear  strain,  the  values  of  G13  reduced  40% 
and  S13  of  75.8  MPa  (11,000  psi)  were  used  to  calculate  £13.  The  resulting 
values  were  2.1%  for  AS4/3501-6  and  1.9%  for  S2  glass/3501-6. 

Substituting  0®  ply  Xj,  (0w  +  0v)»  S13  or  £13  into  eqiiation  8  and  9 

provides  the  stress  concentration  factor  for  each  laminate.  These  values 
are  Usted  in  Table  10,  and  are  seen  to  be  comparable  for  each  laminate 
thickness,  although  they  differ  depending  on  which  compression  failure 
theory  was  used  (see  eq.  8  and  9).  The  fact  that  kt  is  unchanged  with 
thickness  suggests  that  the  strength  of  these  laminates  should  not  change 
with  thickness  and  to  further  investigate  this  observation,  the  values  of  kt 
from  the  48  ply  specimens  were  used  to  predict  the  strength  of  the  96  and 
192  ply  laminates  considering  the  fiber  exit  angle  determined  for  each  of 
these  thicknesses.  These  results  are  shown  in  Figures  7  and  8  and  show  a 
decreasing  trend  in  theoreticed  compression  strength  of  the  same 
magnitude  as  the  experimentally  measured  values.  These  figures  also 
show  that  the  consideration  of  initial  waviness  alone  will  not  change 
compression  strength  with  increasing  specimen  thickness. 

It  shotild  be  noted  here  that  the  stress  concentration  factors  and  the 
distance  over  which  integration  was  performed  to  determine  qn  cannot  be 
applied  in  general  to  the  compressive  failure  of  composite  materials. 
Instead  these  terms  have  been  used  within  in  a  consistent  theoretical 
framework  to  determine  the  effect  of  a  specific  fixture  geometry  on 
compression  coupons  of  increasing  thickness.  For  variations  of  fixture  and 
specimen  geometry  a  comparable  set  of  experiments  with  associated 
analysis  should  be  performed  to  determine  these  factors.  Within  the  scope 
of  this  discussion,  the  above  analysis  has  been  used  to  show  that  a  shear 

U 


based  failure  in  compression  specimens  has  occurred  in  the  region  of  high 
fiber  misalignment,  and  is  related  to  the  shear  stress  caused  by  fiber 
misalignment.  This  correlation  between  theoretical  predictions  and 
experimental  observations  is  otherwise  difficult  to  make  due  to  the 
sensitivity  of  compression  response  to  the  many  local  mhomogeneities  that 
may  be  present. 

When  considering  the  use  of  thick  composites  in  structural 
applications,  the  effect  of  through-thickness  Poisson  expansion  must  be 
closely  considered  in  areas  where  through-thickness  restraint  is  present. 
Through-thickness  restraint  by  test  fixtures  is  beneficial  in  preventing 
premature  end  failiires  in  end-loaded  compression  coupons,  however  the 
effect  of  this  restraint  near  the  gage  region  has  been  shown  to  be 
imdesirable  in  this  study,  and  worsens  as  specimen  thickness  increases. 

In  addition  to  the  compression  test  methods  analyzed  in  this  program,  this 
method  of  accounting  for  the  effect  of  through-thickness  displacements  on 
compression  strength  should  be  extended  to  structural  details  such  as 
joints.  As  with  compression  tests  specimens,  the  tradeoffs  between  desired 
effects  of  restraint  (prevention  of  end  initiated  failures  or  failure  due  to  high 
through-thickness  stresses)  and  imdesirable  effects  must  be  balanced. 

Condusioiis 

This  paper  describes  an  analysis  performed  on  thick  composite 
specimens  tested  to  failure  in  compression.  The  resiilts  have  shovm  that 
there  is  no  inherent  effect  of  thickness  on  the  strength  of  carbon  and  S2 
glass  reinforced  laminates  even  though  a  drop  in  strength  was  observed 
experimentally.  Decreases  in  compression  strength  observed  for 
specimens  from  6.4  to  25.4  mm  have  been  shown  to  be  caused  by  through¬ 
thickness  fixture  induced  effects.  The  fixture  restraint  on  through¬ 
thickness  Poison  expansion  results  in  through-thickness  fiber 
misalignment  that  is  significant  enough  to  account  for  the  reduction  in 
measured  compression  strength. 

This  and  other  results  from  this  investigation  suggest  that  what  has 
been  learned  from  the  evaluation  and  analysis  of  thin  composites  in 
compression  can  be  used  in  the  evaluation  and  analysis  thick  composites  in 
compression.  And  conversely,  what  has  been  learned  from  experiments 
and  analysis  of  thick  composites  can  be  used  to  further  guide  the 
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development  of  failure  theories  for  composite  materials  in  more 
conventional  thicknesses.  In  particular,  the  correlation  of  experimental 
results  with  results  predicted  by  kink  band  based  failtire  theories  suggest 
continued  work  in  this  area.  Current  failure  theories  that  are  based  on  the 
kink-band  failure  mechanism  are  very  simplistic  and  could  benefit  from 
further  development .  The  inclusion  of  nonlinear  shear  stress  strain 
response  in  these  theories  is  an  area  in  need  of  development. 

A  correlation  between  fiber  misalignment  angle  and  compression 
strength  was  shown  both  experimentally  and  theoretically  in  this  work. 
Little  information  exists  on  the  magnitude  and  dispersion  of  misalignment 
in  typical  fiber  reinforced  composites  and  this  area  should  be  further 
investigated.  Data  documenting  the  amplitude  and  wavelength  of  fiber 
waviness  is  almost  nonexistent  and  is  needed  for  further  imderstanding  of 
this  effect  on  compression  strength. 

Finally,  with  the  interest  in  using  composite  materials  for  larger  and 
thicker  structures  the  need  to  account  for  Poisson  induced  fiber 
misalignment  as  documented  in  this  study  becomes  important.  In  many 
joint  concepts  that  could  be  envisioned  for  thick  composite  structures,  a 
through-thickness  restraint  is  present,  and  the  effect  of  this  restraint  must 
be  considered  when  the  material  is  subjected  to  inplane  compressive  stress. 
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Table  1  Efifective  Gage  Section  Expansion. 


Eq.  1 
Free 

Expansion 
mm  (in.) 

Eq.  2 
Preload 
Contraction 
mm  (in.) 

Eq.3 

Restrained 
Expansion 
mm  (in.) 

Eq.4 
Efifective 
Expansion 
mm  (in.) 

AS4/3501-6 
48  ply 

0.0483 

(0.0019) 

0.0076 

(0.0003) 

0.0254 

(0.0010) 

0.0305 

(0.0012) 

96  ply 

0.0838 

(0.0033) 

0.0178 

(0.0007) 

0.0508 

(0.0020) 

0.0508 

(0.0020) 

192  ply 

0.1524 

(0.0060) 

0.0330 

(0.0013) 

0.0965 

(0.0038) 

0.0889 

(0.0035) 

S2/3501-6 

48  ply 

0.0787 

(0.0031) 

0.0025 

(0.0001) 

0.0559 

(0.0022) 

— 

0.0253 

(0.0010) 

96  ply 

0.1372 

(0.0054) 

0.0076 

(0.0003) 

0.1041 

(0.0041) 

1  192  ply 

0.2540 

(0.0100) 

0.0178 

(0.0007) 

0.1981 

(0.0078) 

0.0737 

(0.0029) 
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Table  2  Theoretical  and  Experimental  Laminate  properties  Used  ^n  Finite 
Element  Analysis 


Ex  GPa 
(Msi) 


Ey  GPa 
(Msi) 

Ez  GPa 
(Msi) 


Gxy  GPa 
(Msi) 


Gyz  GPa 
(Msi) 


41.9 

(6.08) 

13.1 

(1.90) 


0.613 


6.00 

(0.87) 


Gxz  GPa  5.03 

(Msi)  (0.73) 


4.34 

(0.63) 


71.8 

(10.41) 

[3.7]l[26]2 


0.072 
[7.3]  [5] 

0.608 
[3.6]  [14] 


S2  glass/3501-6 
[02/90]ns 
Theor.  Exp. 

38.0 
(5.51) 
[7.2]  [31] 

22.9 

(3.32) 

13.4 

(1.95) 

0.156 
[4.6]  [7] 

0.450 
[5.7]  [14] 

0.617 

6.76 

(0.98) 

5.38 

(0.78) 

4.41 

(0.64) 
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Table  3  Lamina  Input  Data  for  Laminate  Plate  Theory  Calculations. 


AS4/3501-6 

S2  glass/3501-6 

EiGPa 

107.1  (15.54) 

48.6  (7.05) 

(Msi) 

[3.5]l 

[4.7] 

E2  GPa 

8.83  (1.28) 

9.72  (1.41) 

(Msi) 

[3.9] 

[3.5] 

E3  GPa 

8.83  (1.28)2 

9.72  (1.41) 

(Msi) 

V12 

0.352 

0.311 

[1.2] 

[5.0] 

Vi3 

0.350 

0.331 

[2.1] 

[5.7] 

V23 

0.649 

0.690 

[2.0] 

[3.4] 

Gi2  GPa 

6.00  (0.87)3 

6.76(0.98) 

(Msi) 

Gi3  GPa 

6.00  (0.87)4 

6.76  (0.98)  I 

(Msi) 

G23  GPa 

3.79  (0.55)5 

3.79(0.55) 

(Msi) 

ai 

m/m/°C 

-1.4  xlO'*^ 

1.3  X  10-5 

02 

m/m/°C 

0.8  X  10-5 

1.1  X  10-5 

as 

m/m/°C 

0.8  X  10-5 

1.1  X  10-5 

1  coefficient  of  variation  (%)  4  G13  assumed  equal  to  G12 

2  E3  assumed  equal  to  E2  5  G23  from  reference  (15) 

3  Gi2  determined  from  [±45]2s 
tension  test 
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Table  4  Number  of  Elements  and  Nodes  in  Finite  Element  Model  for  Each 
Specimen  Thickness. 


Element 

width, 

mm, 

(in.) 

48  ply 
Specimens 

96  ply 
Specimens 

192  ply 
Specimens 

elements 

nodes 

elements 

nodes 

elements 

nodes 

0.7950 

(0.0313) 

24 

97 

80 

285 

320 

1049 

0.3962 

(0.0156) 

96 

337 

320 

1049 

1280 

4019 

0.1981 

(0.0078) 

384 

1249 

1280 

4019 

5120 

15713 

0.0991 

(0.0039) 

1536 

4811 

5120 

15713 

— 

0.0660 

(0.0026) 

3456 

10657 

: 

••• 

•— 
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Table  5  Outer  Ply  Exit  Angle  (0v)  as  a  Function  of  Finite  Element  Mesh  Size. 


1  8v 

(degrees)  Values  for  AS4/3501-6  Specimens 

Element  Width,  mm  (in.) 

0.7950 

(0.0313) 

0.3962 

(0.0156) 

0.1981 

(0.0078) 

0.0991 

(0.0039) 

0.0660 

(0.0026) 

48  ply 

0.648 

0.885 

1.17 

1.52 

1.76 

96  ply 

0.796 

1.06 

1.37 

1.77 

— 

192  ply 

0.914 

1.19 

1.54 

— 

— 

0v  (degrees)  Values  for  S2  glass/3501-6  Specimens 

Element  Width, 

mm  (in.) 

0.7950 

(0.0313) 

0.3962 

(0.0156) 

0.1981 

(0.0078) 

0.0991 

(0.0039) 

0.0660 

(0.0026) 

48  ply 

0.560 

0.768 

1.02 

1.33 

— 

96  ply 

0.662 

0.880 

1.15 

1.49 

— 

192  ply 

0.807 

1.105 

1.36 

— 

— 
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Table  6  Values  of  k  and  m  Used  to  Determine  6v 


Values  of  k  and  m  Used  to  Determine  6v 

AS4/3501-6 

S2  glass/3501-6  I 

k 

m 

k 

m 

48  ply 

0.00173 

0.38 

0.00151 

0.38 

96  ply 

0.00219 

0.38 

0.00185 

0.38 

192  ply 

0.00258 

0.38 

0.00230 

0.38 

Table  7  Outer  Ply  Exit  Angles. 


[02/90188 

[02/90]32s 

AS4/3501-6 

Expansion  Gy,  (®) 

0.789 

0.998 

1.18 

Initial  0w,  (“) 

0.660 

0.660 

0.660 

Total  <t>o,  (®) 

1.45 

1.66 

1.84 

(radians) 

(0.0253) 

(0.0289) 

(0.0321) 

S2  glass/3501-6 

Expansion  0,,  (®) 

0.688 

0.843 

1.05 

Initial  Gw,  (”) 

0.820 

0.820 

0.820 

Total  4>o,  (°) 

1.51 

1.66 

1.87 

(radians) 

(0.0263) 

(0.0290) 

(0.0326) 

Table  8  By  as  a  function  of  Gy  nonlinearity  (degrees). 


Table  10  kt  Determined  From  Compression  Failure  Theories. 


1 

[02/90]8s 

[02/90]i6s 

[02/90]32s 

Average 

1  kt  Argon 

1  AS4/3501-6 

1.92 

2.01 

1.92 

1.95 

S2  glass/3501-6 

2.10 

2.03 

2.10 

2.08 

kt  Budiansky 

AS4/3501-6 

1.04 

1.16 

1.15 

1.12 

S2  glass/3501-6 

1.25 

1.25 

1.35 

1.28 
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THICK-SECTION 
COMPRESSION  TEST  METHOD 


Hardened 
Steel  Plate 


7777777777777777777777Z 


Figure  1  Cross  Sectional  View  of  the  Thick-Section  Compression  Test 
Fixture. 
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Tab  Tab 

Figure  2  Representative  Kink  Band  Geometry  From  a  48  ply  6.4  mm 
AS4/3501-6  Specimen. 
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Effective  Gage  Section  Expansion 
and  Outer  Ply  Exit  Angle 


Fif^ure  3 


f  t  H  ^ 

Bolt  preload  and 
restaining  load 


1/2  effective  gage 
section  expansion 


Undeformed 

specimen 


Schematic  of  Effective  Gage  Section  Expansion  and  Outer  Ply 
Exit  Angle. 
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Finite  Element  Mesh  Discretization 

Coarsest  Mesh  (element  width,  w  =  0.795  mm  [0.0313  in.]) 


Figure  4  Finite  Element  Mesh  Discretization 
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Uniform  Displacement 
Boundary  Condition 


Figure  5  Finite  Element  Model  Geometry  Used  to  Determine  Outer  Ply  Exit 
Angle. 
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Compression  Strength,  MPa 


Compression  Strength  Considering 
Fiber  Misalignment  at  Failure  Stress 
[0/0/90]  AS4/3501-6  Laminates 


Thickness,  mm 

Figure  7  Theoretical  and  Experimental  Comparison  of  Compression 

Strength  vs.  Thickness  For  AS4/3501-6.  kt  Determined  From  48 
Ply  Specimens. 
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Compression  Strength,  MPa 


Compression  Strength  Considering 
Fiber  Misalignment  at  Failure  Stress 
[0/0/90]  S2  Glass/3501-6  Laminates 


Thickness,  mm 

Figure  8  Theoretical  and  Experimental  Comparison  of  Compression 

Strength  vs.  Thickness  For  S2  glass/3501-6,  kt  Determined  From 
48  Ply  Specimens. 
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